Numerical calculation of classical and non-classical electrostatic potentials 
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We present a numerical exercise in which classical and non-classical electrostatic potentials were 
calculated. The non-classical fields take into account effects due to a possible non-zero photon rest 
mass. We show that in the limit of small photon rest mass, both the classical and non-classical 
potential can be found by solving Poisson's equation twice, using the first calculation as a source 
term in the second calculation. Our results support the assumptions in a recent proposal to use ion 
interferometry to search for a non-zero photon rest mass. 



I. INTRODUCTION 

Although classical electromagnetism forbids electro- 
static fields inside empty conducting shells, quantum me- 
chanics suggests that small fields might exist. In the 
spirit of Yukawa's particle-exchange theory of forces Q, 
a modified version of Maxwell's equations was derived to 
account for a possible non-zero rest mass of the photon, 
the exchange particle of the Coulomb force . Because a 
finite photon mass would hmit the range of the Coulomb 
force, these equations violate Gauss's law and make these 
fields possible. 

The experimental search for deviations from 
Coulomb's inverse-square law goes back as early as 
1769 HSU. Although no field has been found at 
the sensitivity level of these experiments, based on the 
predicted sensitivity of the experiments an upper limit 
on the photon rest mass has been determined. The 
most recent tests of Coulomb's law were hmited by the 
sensitivity of the voltage-measurement electronics and 
possible back-action of the measurement process on the 
potentials being measured 0, 0| . 

Progress on these experiments has been slow — the 
limit on the rest mass of a photon from these types of 
experiments has decreased by only a factor of 2.5 in the 
last 35 years. We recently proposed to use ion interfer- 
ometry to iinprove this measurement by several orders of 
magnitude In this experiment, a voltage would be 
applied across a concentric pair of conducting cylindri- 
cal shells. A beam of ions traveling through the inner 
shell would be split and recombined using either physi- 
cal gratings or laser beams. A non-zero electric field in 
the shell would induce a phase shift between the arms of 
the interferometer, resulting in a shift in the interference 
pattern. 

While investigating the feasibility of the experiment, 
we performed several numerical calculations of classical 
and non-classical electrostatic potentials in the proposed 
apparatus. In this paper we discuss the methods and 
results of these calculations. 



II. DEFINITION OF THE PROBLEM 

In the proposed experiment, a non-zero electric field 
will be searched for inside the inner of two concentric 



cylindrical conductors held at different voltages. The in- 
ner conductor's end caps would have several smah holes 
to allow passage of the ions and to ahow laser beams 
and/or wires to enter. For the experiment to work, 
the fringing fields "leaking" through these holes must be 
small compared to the non-classical field under study. 
In addition, the calculations in Q assume that the non- 
classical part of the potential between the interferometer 
gratings is approximately that of an infinite cylinder. To 
verify that these conditions are met, we performed sev- 
eral numerical calculations. 

The setup assumed for these calculations is illustrated 
in Fig. nj The inner conducting shell was assumed to be 
a thin-walled 2.6 m long cylinder, 27 cm in radius. It is 
capped on either end with conducting disks, also 27 cm 
in radius. These disks are 20 cm long to reduce fringing 
fields through the holes. This inner shell is surrounded by 
a second conducting tube with end caps. The outer tube 
is 3.06 m long with an inner radius of 30 cm, giving a 3 cm 
clearance between the inner and outer shells on all sides. 
The inner shell was assumed to be at a voltage V relative 
to the outer shell, which was assumed to be grounded. 
In an actual experiment, V will likely be hundreds of 
kV. But because all of the potentials we calculated scale 
linearly with 1/, we set F = 1 V for our calculations, and 
then scaled the results. 

For simplicity we performed our calculations for a 
cylindrically-symmetric geometry, replacing the holes in 
the inner conductor's end caps with radial slices gener- 
ated by rotating a 1 cm hole about the axis at a radius 
of 25 cm. As such, the fringing fields that we calculate 
will be significantly larger than the actual fields in the 
real apparatus, and the calculation should be considered 
a "worst case" estimate. 



III. EQUATIONS FOR CLASSICAL AND 
NON-CLASSICAL POTENTIALS 



To find the classical fringing-field potential (j)c from a 
given set of static boundary conditions we use Laplace's 
equation. For the non-classical field, we use the counter- 
part to Laplace's equation for massive photons: 
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Figure 1: Layout of conductors. This scale figure gives a cut-away perspective view of the conductors which enforce the 
boundary conditions assumed in our calculations. Arrows point out (a) the inner conducting cylinder, (b) the outer conducting 
cylinder, (c) the disk which caps one end of the inner conductor, and (d) one of the three sets of additional conductors inserted 
into our calculation to explore the effect of mounting objects inside the inner conductor. 



The constant fi in this equation is related to the photon 



rest mass rrij by the relation /i 



yC/h. 



For the non-classical calculation, instead of solving for 
, we solved for the deviation from the classical potential 
p — (j)—(j)c. With this definition, Eq.Q] becomes 



(2) 



Laplace's equation states that V^0c — 0, so the first term 
in the above equation is zero. Also, since fi is known to 
be very small, we expect (j) to be approximately equal to 
(pc- This implies that (pp will be very small compared to 
(pc- As such, we can drop the last term. And since we 
don't know a-priori what /i is equal to, we will normalize 
our equation by defining a new parameter (j)s = (pp/^'^. 
Making this substitution and cancelling out the fi^ in 
both terms we get 



(3) 



Equation is simply Poisson's equation with play- 
ing the part of the charge distribution. As such, our 
simulation need only to be able to solve one equation. 



(4) 



To calculate the classical potential (j)c, we simply replace 
4>x with <pc and insert y — into this equation. To calcu- 
late the non-classical part of the potential, we first calcu- 
late (/)c. Then we then replace cpx with <ps, and insert our 
previously calculated values for (pc as the source term y. 
For an axially-symmetric system, Eq.01can be written as 
a two-dimensional equation in cylindrical coordinates: 
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(5) 



where r and z are the radial and axial coordinates. 

Note that if /i ^ 0, a constant potential is not a so- 
lution to Eq. m and we are not free to arbitrarily define 



the outer conductor to be at V„ 



0. But if this con- 



ductor is Earth grounded, due to the huge capacitance of 
the Earth it is reasonable to assume that Vg will remain 



fairly constant as the voltage on the inner conductor is 
changed. In the proposed experiment the inner conduc- 
tor's voltage relative to the outer conductor would be 
periodically reversed. Rather than measuring the field 
inside the conductor, the difference between the field be- 
fore and after the reversal would be measured. In this 
measurement the the unknown voltage offset Vg will can- 
cel. As such, setting = in our calculations still gives 
us meaningful results. 



IV. METHODOLOGY 

We did our calculations using a finite difference 
method, in which the potential is calculated at points 
on a grid. Derivatives are approximated to second order 
from adjacent points on the grid. Our grid was evenly 
spaced in both dimensions, with points separated by a 
distance a. The potential (j)x at each point on the grid 
will be written as (pij and the source term y at each point 
as yij, where i and j are integers labeling the point. We 
will define i = j ~ at the center of the conducting 
shells, such that the actual coordinates of each grid point 
are equal to = ai and zj = aj. Using these definitions, 
Eq. Isl can be approximated by the equation 



2arj 



y^.3■ (6) 



This can be solved for (j>i j in terms of the known quantity 
yi,j and the value of (px at adjacent grid points: 
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(7) 



We began our simulation with an arbitrary value of 
Lj at each grid point. Then using Gauss-Seidel iteration 
I, this equation was evaluated at each point to produce 
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an updated value of After many iterations, 

eventually converged to the correct values to solve Eq. 
To accelerate convergence, we used the successive over- 
relaxation method 

One group of points which requires attention are the 
points for which i — 0. Because r is always positive 
in cylindrical coordinates, is not defined for nega- 
tive values of i. But there is effectively no difference be- 
tween cylindrical and Cartesian coordinates for the row 
of points along the axis. So for these points we calculated 
derivatives using Cartesian coordinates knowing that the 
points directly below the axis and just into or out of the 
two-dimensional grid should have the same value as the 
point directly above each point on the axis. This gives 
the equation 

0oj = ^ (401J + <l>o,j+i + <Po,j-i) ■ (8) 

In addition to axial symmetry, the conductors have 
mirrored symmetry about their center. This allows us to 
throw away all of the grid points with j < 0, cutting the 
number of grid points in half. Doing this requires us to 
treat the j = points differently, because our grid no 
longer contain values for By symmetry we know 

that — (jii^-j. This allows us to replace 4>i^-i with 
4>i^i in Eq. [3 The point 0o,o is a special point, being 
both a member of the r — Q and z — Q groups of points. 
For this point we use Eq.|Hl but substitute 0o,i for 1^0,-1 • 

After every 100 iterations, the program calculated an 
estimated error at each point by evaluating Eq.[2|at each 
point without changing any values on the grid. We de- 
fined the estimated fractional error to be 

where (f)i^j is the actual value at grid point and <i>ij- 
is the value calculated from Eq.[3 If there were no points 
on the grid with a fractional error larger than 10~*, the 
simulation terminated. 

While experimenting with different conductor geome- 
tries, we significantly decreased the required computation 
time by successively reducing the size of the grid in the 
axial direction. This was done by finding columns which 
had already converged near their final value and then 
using these values as the new boundary conditions for a 
problem involving a smaller grid. To do this the program 
used a variable jmax, initialized to the largest j index 
on the grid, and only updated points with j < jmax. 
After each set of 100 iterations the program stopped and 
calculated the fractional error for each of the points with 
j = jmax. If the fractional error at every point in this 
column was smaller than 10~*, the program would re- 
duce jmax by one, thereby reducing the effective size of 
the grid. The program then repeated this process until it 
found a column which contained at least one point who's 
fractional error was larger than the specified value. Then 
the over-relaxation parameter was re-calculated and the 
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Figure 2: Reducing the grid size. In this process the frac- 
tional error of the outermost column of points (the dark grey 
column) is evaluated. If the fractional error for each of these 
points is less than a set value, the grid size is reduced and the 
next column of points (the medium grey column) is evaluated. 
If the error for each point in this column is smaller than the 
set value, the grid is reduced again and the next column (the 
light grey column) is evaluated. This process continues until a 
column is found which has at least one point with a fractional 
error larger than the specified limit. In this example, the 
right-most white column was the first column found to have a 
point with a fractional error above the set limit, and the grid 
was reduced in the direction of the arrow to the dashed line. 



next set of 100 iterations was performed. This process is 
illustrated in Fig. [2 

The error introduced by stepping in should be negli- 
gible; if a column has converged to within a factor / of 
its final value, the error introduced onto other points by 
"freezing" this column should be of order /. Therefore, 
for a grid with columns, the largest fractional error 
introduced anywhere on the grid by this method should 
be on the order of f\fN~z if the errors are assumed to be 
random, and on the order of /TV^ or smaller otherwise. 

Because we have not done a rigorous theoretical study 
of this method, once we had decided on the final geom- 
etry for our conductors we verified our calculations by 
performing additional computations which did not use 
this "stepping in" method. The result of these calcula- 
tions were identical to those done with the stepping-in 
technique to the seven digits of precision saved at the 
end of the calculations. 



V. CLASSICAL FRINGING-FIELD POTENTIAL 

For the classical fringing-field calculation we are 
mainly interested in how the potential inside the inner 
conductor varies from V , the voltage of the inner con- 
ductor. To keep round-off error from completely mask- 
ing these variations, we made use of the fact that classical 
electromagnetism allows us to arbitrarily add a constant 
potential. So rather than solving Laplace's equation for 
(^c, we instead solved for A(^c = </'c — ^ using the same 

equation but different boundary conditions V ow the 

outer conductor and on the inner conductor. This way 
we found small deviations from zero rather than a finite 
potential. 

The results of this calculation are shown in Fig. Ela). 
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Figure 3: The results of the numerical calculations. All calcu- 
lations were done with a grid spacing of a = 1.67 x 10~^ cm. 
Only the region inside the inner conductor is shown. In (a) a 
semi-log plot of the classical fringing-field potential is shown. 
The peak in the upper right-hand corner is at the location of 
the radial slice in the end cap. In (b) the non-classical po- 
tential (jis is plotted. Frame (c) shows the effect of additional 
conducting objects inside the inner shell. 



We have verified that the errors due to a finite grid are 
negligible by performing this calculation with other grid 
spacings. If the grid spacing a is increased by a factor of 
3, A(/)c at the center of the conductors changes by 9%. If 
a is increased by a factor of 1.5, it changes by only 1.3%. 
We fit these three results to a power law, and found that 
the fit is in good agreement with a fourth data set data 
in which a is increased by a factor of 6. Prom this fit we 
estimate an error on the order of 0.3%. 



To verify the validity of our results we used a series so- 
lution to calculate the field inside of the inner conductor 
assuming that the inner conductor was grounded, and 
that the potential inside of the radial slices in the end 
caps was a fixed constant. This calculation showed that 
the potential at the center of the tube would be 4.5 x 10"'' 
times the potential inside the radial slices, in good agree- 
ment with our numerical calculation. We also found a 
series solution for the field between two grounded cylin- 
ders with a fixed potential applied to the gap at the ends. 
The length of the cylinders was assumed to be equal to 
the sum of the lengths of the two end caps, with a gap 
between them equal in size to the radial slice in the end 
caps. The field in the center of the gap turned out to be 
2 X 10^'' times smaller than the potential at the end, also 
in agreement with our numerical calculation. 



VI. NON-CLASSICAL POTENTIAL 

To calculate 0^ , we set y — (l)c in Eq. El and used the 
same methods discussed above. Since the potential on 
the conducting surfaces is given, the deviation from the 
classical potential should be zero on these surfaces. So 
the boundary conditions for the (j)s calculation were that 
0s = on both of the conductors. 

The source term y in this calculation is equal to (j>c. We 
can obtain this by adding V to our already completed cal- 
culation of A(j)c. Inside the inner conductor, where A(j)c 
is very small, this results in values accurate to the total 
precision allowed by the double-precision fioating-point 
format. But between the two conductors, adding V to 
A(j)c results in lost precision. We are only interested in 
the fields inside the inner conductor, which should not be 
affected by this lost precision — the non-classical poten- 
tial is predominantly generated by the local source term 
y rather than from fringing fields generated outside the 
inner conductor. But to be extra careful we also calcu- 
lated 4>c directly using the correct boundary conditions 
(the inner conductor at V and the outer one at 0), and 
used these values wherever 0c was less than V/2. 

The results of this calculation are shown in Fig. E^b) . 
Figure EJc) shows the results of a similar calculation in 
which additional conductors were added to simulate op- 
tics and other objects inside the inner shell. These con- 
ductors consisted of a ring and a cylinder at the three 
locations where gratings would be positioned in the in- 
terferometer. The rings were assumed to be 6 cm wide 
and 1.5 cm thick, just touching the inner shell, and the 
cylinders to be 6 cm long and 1.5 cm in radius, centered 
on the axis of the inner shell, as shown in Fig. These 
conductors were assumed to be at the same potential as 
the inner shell. 

We estimate the error in these calculation to be very 
small — to the seven digits output by our program there 
was no change in 0^ when we increased the grid spacing 
by a factor of 3. 



5 




Figure 4: Comparing the different potentials. The potential 
is plotted as a function of radius at four different points along 
the axis. In (a) the non-classical field is shown. In (b) the 
non-classical field is shown for the case in which additional 
conductors have been placed inside the inner shell. In (c) 
the classical field is shown. The lines marked with circles 
and diamonds are at axial distances of and 100 cm from 
the center of the conductors — the locations of the gratings. 
The lines marked with + and x are at axial distances of 33 
and 67 cm, respectively. All potentials were plotted assuming 
that V = 200 kV and ^J. = 2A x 10"" m"\ The thick grey 
lines in (a) and (b) represent the analytical solution for the 
non-classical field inside of an infinite cylinder. 



VII. IMPLICATIONS 

The fields are compared in Fig. 21 While this figure 
shows the potential at each point moving radially from 



the axis, for the proposed ion interferometer experiment, 
all that is important is the field at the location of the ion 
beam. As such, the most important information in this 
figure is the radial slope of the potential at r = 25 cm. 

In Fig. EJa) we see that the non-classical field is ap- 
proximated very well by the field of an infinite cylinder. 
In Fig. ^b) we see that although the effect of small ob- 
jects inside the shell on the non-classical potential is not 
negligible, it should not greatly change the sensitivity of 
the experiment as long as care is taken. Because ax- 
ial symmetry is assumed, the additional conductors took 
the form of large rings rather than small rectangles which 
would better approximate an optical mount, and this plot 
can be considered an extreme "worst-case" estimate. 

Note that the vertical axis in Fig.^c) is about 2.4 x 10* 
times smaller than in (a) and (b), indicating that the 
fringing-fields inside the inner conductor should be com- 
pletely negligible for values of /i much smaller than 
2.8 X 10~^^ m~^ (corresponding to a photon rest mass 
of 1 X lO^'^^ kg, over 600 times smaller than the cur- 
rent experimental limit measured in Q). Consequently, 
fringing fields should not be a problem in the proposed 
experiment. 



VIII. CONCLUSIONS 

In conclusion, we have conducted a numerical study of 
classical and non-classical electrostatic potentials in an 
axially-symmetric nested conductor configuration. The 
results show that the assumptions in our recently pro- 
posed ion-interferometry experiment are valid. The cal- 
culations show that for values of ^ much smaller than 
the current experimental limit, non-classical fields should 
still dominate over fringing fields. Furthermore, we have 
shown that the non-classical field approximates the sim- 
ple field of an infinitely long set of conductors. 

We would like to acknowledge Ross Spencer for his help 
on every aspect of this study. This work was funded by 
BYU's Office of Research and Creative Activities. 
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